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Water adsorption and dissociation on SrTiO3(001) revisited: 
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We present a comprehensive density-functional theory study addressing the adsorption, dissoci- 
ation and successive diffusion of water molecules on the two regular terminations of SrTiO3(001). 
Combining the obtained supercell-geometry converged energetics within a first-principles thermo- 
dynamics framework we are able to reproduce the experimentally observed hydroxilation of the 
SrO-termination already at lowest background humidity, whereas the Ti02-termination stays free 
of water molecules in the regime of low water partial pressures. This different behavior is traced back 
to the effortless formation of energetically very favorable hydroxyl-pairs on the prior termination. 
Contrary to the prevalent understanding our calculations indicate that at low coverages also the 
less water-affine Ti02-termination can readily decompose water, with the often described molecular 
state only stabilized towards higher coverages. 

PACS numbers: 68.43.Fg 71.15.Mb 68.47.Gh 
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I. INTRODUCTION 

In many of its applications e.g. as photocatalyst, gas 
sensor or growth template SrTiOa is deliberately or un- 
intentionally exposed to aqueous environments. This 
has motivated a large number of studies on fundamen- 
tal aspects of the interaction of water with this partic- 
ular material^. As for many other SrTiOa properties, 
they have most of all revealed a sensitive dependence on 
the detailed surface morphology and therewith prepara- 
tion procedure. Within the Surface Science philosophy 
this shifts the focus to the nominally "ideal" two non- 
polar terminations of SrTiO3(001), cf. Fig. [H aiming 
to establish firm answers for these well defined refer- 
ences. The picture emerging from corresponding stud- 
ies is an enhanced reactivity of water with the SrO- 
termination^^^ . While this termination seems to hydrox- 
ylate readily, weak molecular adsorption appears favored 
at the Ti02-termination in humid environments. Still, 
primarily due to the difficulties of preparing a perfect 
SrTiO3(001) surface that exhibits only one defect-free 
termination, a discomforting degree of uncertainty con- 
cerning the adsorption state, structure and energetics re- 
mains. A key experiment in this respect is therefore the 
Friction Force Microscopy (FFM) study by Iwahori et ai, 
which exploits the FFM scanning capability to explore 
the adsorption of water separately on Ti02- and SrO- 
terminated domains^. For a range of water exposures 
a change of the friction force was only measured above 
SrO-terminated domains and attributed to surface hy- 
droxylation. This confirms the stronger affinity of water 
to this domain and sets an upper bound for the water 
bond strength at the Ti02-termination. 

Several theoretical works have already been performed 
to complement these experimental insights. Wang et al. 
carried out first density-functional theory calculations 
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FIG. 1: (Color online) Schematic illustration of the SrTiOs 
perovskite bulk structure and the two regular terminations 
of the (OOl)-surface. Big red (dark) spheres: strontium, big 
green (bright) spheres: titanium, small blue (bright) spheres: 
oxygen. 



(DFT) with a gradient corrected exchange-correlation 
(xc) functional and determined a binding energy of 
molecular water at the Ti02-termination in somewhat 
the range indicated by experiment^. Unfortunately, the 
calculations were restricted to small surface unit-cells 
(thereby modeling a high overlayer density), the Ti02- 
termination, and did not compare dissociative and molec- 
ular adsorption modes. Using a hybrid xc functional 
Evarestov, Bandura and Alexandrov expanded on the lat- 
ter two points^. However, in their calculations molecular 
and dissociative adsorption mode were essentially degen- 
erate at the SrO-termination and only marginally more 
stable than molecular water at the Ti02-termination - 
thereby questioning the prevalent understanding of a 
more facile surface hydroxylation of the prior termina- 
tion. Very recently, Baniecki et al. showed that the sim- 
ilar stability at the two terminations is an artefact of the 
employed small surface unit-cells^. However, their calcu- 
lations in larger unit-cells did not systematically compare 



dissociated and molecular adsorption states. 

This situation motivates us to revisit the problem with 
a comprehensive DFT study that systematically adresses 
the adsorption, dissociation and decomposition of water 
molecules at both regular SrTiO3(001) terminations. We 
show that apart from the coverage-dependence already 
identified by Baniecki et at, also the employed slab thick- 
ness is a hitherto not sufficiently appreciated factor for 
the binding energetics. The finally obtained converged 
results are in full agreement with the view deduced from 
experiment, albeit with a small additional twist: Both 
terminations are able to dissociate water at low coverage 
without appreciable barrier. The molecular adsorption 
state discussed at the Ti02-termination is only stabilized 
at higher coverages. At the SrO-termination the dissoci- 
ated water stabilizes in hydroxyl pairs, which we identify 
as the main factor behind the higher water affinity of this 
termination observed in the FFM measurements. 
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FIG. 2: (Color online) Perspective view of the most favorable 
adsorption geometry Al at the SrO-termination, in which the 
water dissociates into an adjacent hydroxyl pair. Atoms dis- 
cussed in the text and Table [T] are labeled. Coloring here and 
in all consecutive figures follows the one of Fig. [TJ with small 
black spheres denoting hydrogen atoms. 



All DFT calculations were performed with the 
CASTEP^" code using a plane- wave basis together with 
ultrasoft pseudopotentialsii as provided in the default li- 
brary, and the GGA-PBE functional^^ to treat electronic 
exchange and correlation. Adsorption at the two regu- 
lar SrTiO3(001) terminations shown in Fig. [T]was mod- 
eled in supercell geometries using mirror-symmetric slabs 
based on the GGA-PBE optimized bulk lattice constant 
(flo = 3.938 A), with adsorption at both sides and a vac- 
uum separation exceeding 13 A. As further detailed be- 
low long-range geometric relaxation effects lead to a slow 
convergence of the binding energetics with slab thick- 
ness. We therefore employed up to 11 layer slabs, in 
which the central three atomic layers were always fixed 
at their respective bulk positions in the cubic (Pm3m) 
perovskite phase, while all other substrate atoms were 
fully relaxed until the absolute value of all corresponding 
forces dropped below 0.05 eV/A. 

The central energetic quantity taken from the DFT 
calculations is the binding energy with respect to gas- 
phase water, defined as 



to the experimental 0.96 A and — 2.52eV, respectivelyi^. 
Systematic tests in (1 x 1) surface unit-cells indicate that 
a cut-off energy for the plane- wave basis set of 430 eV and 
(8x8x1) Monkhorst-Pack (MP) grids^^ for the Brillouin 
zone integrations ensure a numerical convergence of 
within ±20meV. For calculations in larger surface unit- 
cells the MP grids were reduced to maintain the same 
sampling of reciprocal space. 

Energetic barriers and minimum energy paths (MEPs) 
have been computed using the nudged elastic band 
(NEB) methodic in the implementation of the "Atomic 
Simulation Environment""'^''. Corresponding calculations 
relied on seven atomic layer slabs and (3 x 3) surface 
unit-cells throughout. Systematic tests confirm that nei- 
ther for the obtained transition nor initial and final state 
geometries spin-polarization plays a role. Furthermore, 
in all conformations relevant for this work no artificially 
delocalized electronic states were encountered that could 
result from the employed semi-local xc functional and 
which have for example been clearly demonstrated for 
the case of hydrogen impurities in bulk Ti02»^ 
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Here i?H2 0@surf is the total energy of the adsorbate cov- 
ered surface (either with molecular or dissociated water), 
-Bsurf the total energy of the clean surface, and £'H2 0(gas) 
the total energy of the gas-phase H2O molecule (all three 
computed at the same plane- wave cutoff) . The factor 1/2 
accounts for the fact that adsorption is at both sides of 
the slab and in the sign convention of eq. ([T]) a negative 
value of the binding energy indicates that adsorption is 
exothermic. The reference energy of the isolated water 
molecule i?H20(gas) was calculated in a (20 X 20 X 20) A 
super-cell resulting in an optimized OH-bond length of 
0.99 A and a dissociation energy of —2.47 eV as compared 
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Testing all high-symmetry sites offered by the regular 
SrO-termination we obtain a clear energetic preference 
for the adsorption geometry Al displayed in Fig. [2] and 
further quantified in Table HI In this geometry the water 
molecule dissociates into two adjacent hydroxyl groups. 
The 01 and HI atom from the dissociated water gen- 
erate a protruding hydroxyl group. Here, the 01 atom 
sits essentially in a bridge site between two surface stron- 
tium cations Sri and Sr2, which is the position it would 
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Al 




Bl 


B2 


Hl-Ol 


0.97 A 


Hl-Ol 


0.98 A 


0.98 A 


H2-02 


1.01 A 


H2-02 


1.73 A 


0.98 A 


H2-01 


1.60 A 


H2-01 


1.01 A 


2.78 A 


Ol-Srl 
01-Sr2 


2.59 A 
2.55 A 


Ol-Til 


2.21 A 


1.84 A 


02-Til 


2.28 A 


Til-03 


1.94 A 


2.29 A 



TABLE I: Selected bond distances in the dissociated adsorp- 
tion geometry Al at the SrO-termination, and the molecular 
Bl and dissociated B2 state at the Ti02-termination. The 
labels for the individual atoms are defined in Figs. [2]and[Sl 
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A'alab 


(ixl) 


(2 X 2) 


(3 X 3) 




Al 


Al 


Al 


7 


-0.83 


-1.18 


-1.22 


9 


-0.84 


-1.24 


-1.30 


11 


-0.83 


-1.28 


-1.35 


5 


-0.79^ 






7 


-0.811 


-1.10=^ 





TABLE II: Binding energies of one water molecule per sur- 
face unit-cell in the most favorable adsorption state Al at the 
SrO-termination shown in Fig. [2] and as a function of the 
number of slab layers A^siab employed. Also shown are the 
values computed by Evarestov et al& using a hybrid xc func- 
tional and by Baniecki et al. using a gradient-corrected xc 
functional similar to the one employed here. All values in eV. 

also take in a continuation of the perovskitc lattice. The 
split-off H2 atom of the adsorbing water forms a second 
hydroxyl group together with a lattice oxygen anion 02. 
The two hydroxyl groups are strongly tilted towards each 
other, suggesting the formation of a hydrogen-bond that 
we will further qualify below. Corresponding bonds with 
a length comparable to the (ioi-H2 — 1.6 A determined 
here have recently also been reported for hydroxyl ad- 
sorption geometries on alkaline-oxide surfaces^®. 

The hydroxylation lifts the lattice 02 anion quite 
strongly out of the surface plane. This increases the bond 
length to the underlying Til atom by 19%. In turn, the 
latter reinforces its bond to the underlying third layer O 
anion with a reduction of the corresponding bond length 
by 25%. This long range geometric relaxation sequence 
is the major cause for a slow convergence of the binding 
energy with the number of slab layers employed in the 
calculations. As detailed in Table HIl even for 11 layers 
slabs there are still remarkable changes of of the order 
of 50 meV in the larger surface unit-cells. As the geomet- 
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FIG. 3: (Color online) Binding energy, internal 01-H2 bond 
length and adsorption induced change of the Hirshfeld charges 
as a function of the vertical height z of the center of mass of 
an approaching water molecule above the SrO-termination. 
The different atom labels follow the definition given in Fig. 
[2] The zero reference for the vertical height corresponds to 
the equilibrium height in the adsorbed state Al. 



ric relaxations can not properly develop in the smallest 
(1 X 1) cells, there is no such dependence there. On the 
one hand this reveals how misleading slab-thickness con- 
vergence tests in such cells can be, on the other hand 
it allows us to compare our results to those reported in 
the earlier studies by Evarestov et al^ and Baniecki et 
al.^ , which employed corresponding thin slabs in their 
calculations restricted to smaller cells. As apparent from 
Table [IT] (and with equivalent findings compiled in Table 
mil for the Ti02-termination) the numbers compare all 
very well, despite the use of a hybrid xc functional in the 
study by Evarestov et al^. We take this as confirmation 
that also the gradient-corrected xc functional employed 
here is able to capture the essential physics - at a (still) 
significantly lower computational cost that in turn en- 
ables us to perform calculations in much larger surface 
unit-cells. The data for these cells also compiled in Ta- 
ble |TT] on the other hand demonstrates that the latter is 
a crucial point as the binding energy exhibits a strong 
coverage dependence, reflecting overall repulsive interac- 
tions consistent with the TPD data from Wang et al^. 
A value fairly representing the low-coverage limit is only 
reached in (2 x 2) cells, and as we will discuss in Section 
III.C below, it is this limit that is the appropriate one to 
discuss the Iwahori FFM experiments^. 

At these relevant low coverages the dissociated water 
state Al can be reached without any energetic barriers. 
Figure [3] shows the evolution of the binding energy as 
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a function of the molecule's vertical height above the 
surface. The barrier less increase of the binding energy 
upon approaching the surface is accompanied by a grad- 
ual charge transfer to the molecule as reflected by the 
computed Hirshfeld charges^^ also displayed in Fig. [S] 
The adsorption process fits therefore perfectly into the 
picture for dissociative adsorption of water, with the wa- 
ter as a Lewis acid and the surface, primarily the 02 
anion, acting as a Lewis baseiii^. 

While the dissociation into the adjacent hydroxyl pair 
is thus non-activated, this is distinctly different for the 
consecutive decomposition process. Judging from the 
stability of each hydroxyl group alone, the initially re- 
quired breaking of the hydroxyl pair will occur via a 
diffusion hop of the protruding 01-Hl group. Figure |3] 
shows the computed energy profile for this disintegration 
pathway. Starting from the most favorable geometry Al 
a sizable barrier of Ai?Ai2 = 0.93 eV needs to be sur- 
mounted. This leads to a shallow metastable geometry 
A2, in which the 01-Hl group sits in a perovskite-type 
bridge site one unit-cell further away than before. With 
a computed barrier in excess of 1 eV for consecutive diffu- 
sion of the surface H2 atom, the most favorable pathway 
for further decomposition proceeds instead through the 
reorientation of the 01-Hl hydroxyl- group to point away 
from the original hydroxyl partner. This flip is only con- 
nected with a negligibly small barrier A£'a23 = 0.08 eV, 
indicating a very fast dangling dynamics. The ensuing 
diffusion of the separated 01-Hl hydroxyl-group through 
hops to adjacent binding sites is then characterized by an 
intermediate barrier A£^a34 — 0.41 eV. Here, we notice 
that the energetic minima corresponding to the equiv- 
alent states A3 and A4 are aligned on the same level. 
This suggests no further significant long-range interac- 
tion between the immobile 02-H2 group and the moving 
01-Hl hydroxyl. On the other hand, this energetic level 
is located 0.67 eV above the one of the initial configura- 
tion Al, which thus directly reflects the additional bond 
strength resulting from the hydroxyl pairing. With such 
strong attraction between the two hydroxyls generated 
from the dissociation of one water molecule, the percent- 
age of freely diffusing hydroxyl groups or H atoms result- 
ing from this process will be very small up to very high 
temperatures. 



B. Ti02-termination 

In agreement with the preceding theoretical studies^"^ 
we also determine a molecular adsorption state on the 
Ti02-termination. The adsorption geometry of this state 
henceforth denoted Bl is depicted in Fig. [SJ with se- 
lected bond distances compiled in Table HI Adsorption 
in this mode induces only small geometric relaxations of 
the SrTiO3(001) substrate, and accordingly we only ob- 
serve a weak dependence of the computed binding en- 
ergy with the number of layers employed in the slab 
model, cf. Table Hill The coverage dependence is equally 
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(ix 
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(2 X 


2) 


(3 X 


3) 




Bl 


B2 


Bl 


B2 


Bl 


B2 


7 


-0.77 


-0.58 


-0.74 


-0.85 


-0.71 


-0.90 


9 


-0.78 


-0.58 


-0.73 


-0.92 


-0.72 


-0.97 


11 


-0.78 


-0.59 


-0.73 


-0.95 


-0.72 


-0.99 


5 


-0.73^ 


-0.64^ 










7 


-0.73a 




-0.74a 









TABLE III: Binding energies of one water molecule per sur- 
face unit-cell in the two most favorable adsorption states at 
the Ti02-termination shown in Fig. [5] and as a function of 
the number of slab layers A'siab employed. Also shown are the 
values computed by Evarestov et aZ.— using a hybrid xc func- 
tional and by Baniecki et al. using a semi-local xc functional 
comparable to the one employed here. All values in eV. 



weak and in contrast to the dissociated state at the SrO- 
termination slightly attractive. We also compute a non- 
activated adsorption pathway into this molecular bound 
state, which is essentially characterized by an increas- 
ing charge depletion around the H2 atom to stabilize the 
01-Til molecule-surface bond. 

As illustrated in Fig. |6] the full separation of the H2 
atom to reach a dissociated state B2 is accompanied by 
a small energy barrier Ai<^Bi2 = 0.09 eV. In this re- 
spect, the molecular state Bl can be considered as a 
true precursor and its stabilization (in contrast to the 
SrO-termination) flts well to the expectation of a weaker 
acidity of the Ti cations compared to the Sr cations on 
SrTiO3(100)^".. On the other hand, at the low coverage 
corresponding to the (3 x 3) surface unit-cell calculation 
behind Fig. [6] the molecular state Bl exhibits only a 
weakly pronounced metastability, with the ensuing disso- 
ciated state B2 (further characterized in Fig. [5]and Table 
H]) more favorable by more than 0.2 eV. Intriguingly and 
as detailed in Table Hill this energetic ordering reverses at 
higher coverages, with the molecular state Bl more stable 
in the smallest (1x1) surface unit-cell. We therefore fully 
reproduce the preference for molecularly adsorbed water 
obtained in the previous calculations-"", but show that 
this preference is restricted to rather dense packings. In 
fact, at these coverages one may well imagine an even fur- 
ther increased stability of mixed phases containing both 
molecular and dissociated water molecules compared to 
the two pure phases studied here. We did not pursue cal- 
culations along this line though, as we will demonstrate 
in the following Section that it is the regime of lower cov- 
erages that is e.g. relevant for the Iwahori FFM experi- 
ments, and in this low coverage regime dissociation into 
two surface hydroxyl groups is the preferred adsorption 
mode with only a slight activation barrier to overcome. 
We note that these flndings are in no contradiction to the 
frequent interpretation of experimental data exclusively 
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FIG. 4: (Color online) Energy profile for surface diffusion of the protruding Ol-Hl hydroxyl group over the SrO-terminated 
surface. The initial state Al corresponds to the hydroxyl-pair after the immediate water dissociation process shown in Fig. [S] 
The transition from state Al to A2 breaks the hydroxyl-pair and is accompanied by a substantial energetic barrier. Config- 
urations A2 and A3 differ essentially in the mutual orientation of the two hydroxyl-groups, whereas in state A4 the diffusing 
species has further increased the distance to its original position Al by hopping to the adjacent binding site. 




FIG. 5: (Color online) Perspective view of the molecular Bl (left) and dissociated B2 (right) adsorption geometry at the 
Ti02-termination. Atoms discussed in the text and Table |T] are labeled. In a second dissociated state (not shown), which is 
energetically equivalent to the state B2, the Ol-Hl group is rotated by 180° degrees and thus points in the same direction as 
the tilted 02-H2 hydroxyl. 



in terms of molecularly adsorbed water—""—. Even though 
an absolute coverage calibration was rarely achieved, it 
is quite clear that these studies operated mostly in the 
higher coverage regime. Moreover and as already dis- 
cussed by Baniecki et alM- frequently employed finger- 
prints for "molecular" water (like the appearance of two 
adsorbate peaks below the 02p valence band in ultravi- 
olet photoemission spectroscopy (UPS)"*'^) are not nec- 
essarily unambiguous. Indeed, we have verified that the 
valence states due to the hydroxyl groups in the B2 ad- 
sorption mode lie at essentially the same energy as the 
3ai molecular orbital of the adsorbed water in the Bl 



mode, such that even a larger fraction of dissociated wa- 
ter would not easily be distinguished in the UPS spectra. 

In contrast to the SrO-termination the preferred 
mechanism for a further disintegration of the adjacent 
hydroxyl-groups of the dissociated state B2 is not via 
hopping of the protruding Ol-Hl group. We compute 
a large barrier of 1.5 eV for a corresponding hop to a 
neighboring Ti cation. Instead, continued diffusion is at 
this termination much more effective via surface hopping 
of the split-off 112 atom. As illustrated in Fig. [7] this 
mechanism proceeds through a sequence of hops with a 
barrier of 0.51 eV and reorientation flips with a barrier of 
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FIG. 6: (Color online) Binding energy, internal 01-H2 bond 
length and change of the Hirshfeld charges along the reaction 
pathway from the molecular bound state Bl to the dissociated 
state B2 at the Ti02-termination. The different atom labels 
follow the definition given in Fig. [2] The zero reference for 
the Hirshfeld charges corresponds to gas-phase water and the 
clean surface. 



0.38 eV. In the hops the H2 atom breaks the bond to its 
directly coordinated lattice O partner and jumps to a di- 
rectly adjacent lattice O anion. In the consecutive flip the 
thereby formed hydroxyl group changes its tilt direction 
to the opposite side, therevifith enabling another hop to 
lattice O anions now belonging to a neighboring surface 
unit-cell. Interestingly, the rate-limiting step of 0.51 eV 
for this diffusion mechanism agrees almost perfectly to 
the calculated activation barrier for proton diffusion in 
bulk SrTiOjr^. A disintegration of the hydroxyls gener- 
ated through the dissociative adsorption of water is thus 
not as limiting as at the Ti02-termination and might 
actually proceed both through surface and sub-surface 
proton diffusion. 



C. First-principles thermodynamics 

In order to directly address the Iwahori FFM experi- 
ments we now combine the supercell-converged energetics 
for adsorbed water at the two terminations and at differ- 



ent coverages within a first-principles atomistic thermo- 
dynamics framework22ri2£. In this approach we assume 
the surface to be in equilibrium with a surrounding water 
vapor environment characterized by a chemical potential 
/^H2 0- For each surface termination the stable surface 
structure at a given /iH2 then minimizes the surface free 
energy, defined as 

7(/^H20) = [GHaOQSurf — Gsurf — ^H20MH2o] • (2) 

Here, GH20@Surf is the Gibbs free energy of the surface 
covered with iVnaO water molecules per surface unit-cell, 
Gsurf is the Gibbs free energy of the corresponding clean 
surface, and 7(/iH2o) is normalized to energy per unit 
area by dividing through the surface area A of the em- 
ployed surface unit-cell. Aiming only for a first assess- 
ment we approximate the difference (GnaOesurf — Gsmf) 
with the difference of the corresponding total energies, 
thereby neglecting vibrational free energy and configu- 
rational entropy contributions. As discussed in detail in 
ref. (|26[ ) this is largely justified, as most of these con- 
tributions effectively cancel in the difference. A notable 
exception to this arises from the vibrational free con- 
tributions from the adsorbate functional groups, which 
particularly for hydroxyl or water groups are in general 
not negligible^^. We nevertheless omit them here, real- 
izing from the results presented below that the relevant 
surface terminations to discuss the Iwahori FFM experi- 
ments correspond to rather low adsorbate densities. 

Contact to the experimental environments can be 
made by exploiting the relation between chemical poten- 
tial and gas-phase temperature and pressure. For this we 
first separate the total energy contribution to the chem- 
ical potential 

A/XH2o(r,PH2o) = MH2o(7',PH2o) - £^H20(gas) • 

At standard pressure p° = 1 bar the relative 
A/iH2o(2^5PH2o) then be derived from enthalpy H 
and entropy S differences tabulated in thermochemical 
tables2^ 

AAiH2o(T,PH2o) = (^(7',Ph2o) - ■H^(0K,Ph2o) " 

T {S{T,p°^.^o)-S{OK,p°^^o)) ■ 

From this the relative chemical potential at any other 
pressure follows finally from the ideal-gas relation 

A^H2 0(I',PH2 0) = A^H20(r,PH20) + ^BT ( ) , 

where fee is the Boltzmann constant. The latter rela- 
tion is appropriate as long as the discussion is restricted 
to humid environments, i.e. those corresponding to the 
gaseous state of water for any temperature or pressure. 
This sets an upper bound for A^h^q — —0.91 eV, which 
corresponds to the chemical potential of water at the (ex- 
perimental) critical point^. As we will see in the follow- 
ing the experimental conditions employed in the Iwahori 
FFM experiments all fall below this Il2 0-rich limit. 
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FIG. 7: (Color online) Energy profile for surface diffusion of the split-off H2 atom over the Ti02-terminated surface. The initial 
states Bl and B2 correspond to the molecular precursor and dissociated hydroxyl-pair shown in Fig. [5l Further disintegration 
takes place via a hop of the surface H2 atom to an adjacent lattice O anion (B3) with a barrier AiJB23 = 0.51 eV, followed 
by a flip of the orientation of the newly formed hydroxyl group towards the neighboring surface unit-cell (B4) with a barrier 
Ai5B34 ~ 0.38 eV. The ensuing equivalent hop to an even more distant lattice O anion (B5) exhibits essentially the same energy 
profile, indicating overall only small longer-range lateral interactions between the hydroxyl groups. 
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FIG. 8: (Color online) Surface free energies, cf. eq. ((2|, of the most stable water adsorption geometries at the different coverages 
at the SrO-termination (left) and Ti02-termination (right). In the top x-axis, the dependence on the water chemical potential 
is converted into pressure scales at 300 K and 500 K. The plain (blue) background boxes mark the region of gas-phase conditions 
probed by Iwahori et al. in their FFM experiments^, while the dotted (gray) boxes indicate the region above the H20-rich 
limit, i.e. where the present approach assuming equilibrium with water vapor is no longer strictly applicable. 



Figure [8] summarizes the obtained results for the two 
terminations. Starting the discussion with the SrO- 
termination, the graph shows the surface free energies 
obtained for the dissociated hydroxyl-pair state Al at 
the three coverages corresponding to one water molecule 
per (3 X 3), (2 x 2) and (1 x 1) surface unit-cell. Reflect- 
ing the absence of significant longer-range lateral inter- 



actions between different hydroxyl-pairs at the two lower 
coverages, the corresponding structures become more sta- 
ble than the clean surface within a narrow chemical po- 
tential range. As apparent from Fig. [S] the experi- 
mental gas-phase conditions probed in the FFM experi- 
ments by Iwahori et al.^ (10^^^ bar < pn^O < 10~^bar, 
T — 300 K) fall clearly above this threshold. The calcu- 
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lations therefore fully support the authors' surmise that 
the observed increase of the friction-force coefficient on 
the SrO-terminated domains is due to surface hydroxy- 
lation. The observation that the Ti02-termination stays 
free of water molecules in the same probed environments 
is also fully consistent with our calculations. The much 
smaller binding energy for the dissociated state B2 that 
is most stable at the lower coverages yields a crossing 
of the corresponding (3 x 3) and (2 x 2) lines with the 
clean surface reference in Fig. [5] only at chemical poten- 
tials towards the upper end of the experimental range. 
While our calculations thus confirm the interpretation of 
the FFM data, the consistency demonstrates vice versa 
that the supercell-converged absolute binding energetics 
from the employed semi-local xc functional, which stands 
behind the thermodynamics lines in Fig. [SJ is quite reli- 
able. This concurs with the earlier observation that the 
semi-local energetics agrees closely to the one obtained 
by Evarestov et al. with a hybrid xc functional and re- 
confirms our assessment that the discrepancies between 
experiment and theory existing before this work were less 
due to the treatment of electronic xc, but more due to 
the use of restricted supercell geometries. 

At both terminations the dense (1x1) overlayers that 
were at the focus of the preceding theoretical works be- 
come most stable only at a very high chemical potential 
of A/iHjO ~ — 0.7eV. This is far outside the relevant 
conditions for the Iwahori experiments and in fact above 
the H20-rich limit, where the present thermodynamic 
approach is no longer strictly valid. This said it is never- 
theless intriguing to note that recent FFM experiments 
by Kato et al. employed much higher water pressures 
in excess of 10"^ bar at room temperature^^.. They ini- 
tially observed a loss of the Iwahori friction-force con- 
trast between the two SrTiO3(001) terminations, which 
reappeared at a pressure of 10~^ bar. The latter thresh- 
old corresponds to a A/iHaO = — 0.6eV, which is only 
by about 0.1 eV smaller than where our calculations 
would predict the stabilization of the dense (1x1) over- 
layers, and specifically the hydroxylated Al phase on 
the SrO-termination and the molecular Bl phase on the 
Ti02-termination. In contrast, for lower chemical poten- 
tials down to A^HaO = — 0.77eV (corresponding to the 
lower pressure limit of 10^^ bar employed by Kato et al.) 
both terminations would simply be hydroxylated with 
less dense arrangements. If hydroxylation is indeed the 
contrast mechanism behind the FFM observations, this 
would be fully consistent with the interpretation in terms 
of the formation of a condensed water layer at the thresh- 
old pressure of 10^^ bar at room temperature proposed 
by Kato et al^. 

IV. SUMMARY AND CONCLUSIONS 

In summary our comprehensive DFT calcula- 
tions demonstrate that both regular terminations of 



SrTiO3(100) are initially able to dissociate water without 
appreciable barrier. At the SrO-termination this leads 
to stable hydroxyl-pairs in the entire sub-monolayer 
regime. At the Ti02-termination a molecular adsorption 
state is a metastable precursor at low coverages and 
becomes the most favorable adsorption mode in a dense 
overlayer. Combining the computed supercell-converged 
energetics for these adsorption structures within a 
first-principles atomistic thermodynamics framework 
we can fully rationalize the friction-force microscopy 
experiments by Iwahori et al., which indicate a much 
higher water affinity of the SrO-termination. The 
absolute binding energetics obtained with the semi-local 
exchange-correlation functional employed in our study 
is therefore fully consistent with the stringent bounds 
set by these measurements. The detailed analysis of 
the continued disintegration of the dissociated water 
molecules reveals a strong pairing mechanism of the two 
surface hydroxyls generated from the dissociation of one 
water molecule at the SrO-termination. Rather than the 
stronger acidity of the Sr cations in SrTiO3(001) it is 
the additional stabilization due to this hydroxyl pairing 
which stands behind the notably different water affinity 
of the two terminations. A similar pairing mechanism 
has recently been reported for the alkaline-oxides CaO 
and BaO^^ and we believe this cooperative effect to 
be an important general feature for low-hydrated oxide 
surfaces. Another important observation is that the O 
atom of the protruding hydroxyl group that has formed 
as a result of the water dissociation process sits in 
both terminations at the site it would also take in a 
continuation of the perovskite lattice structure. This 
is distinctly different to adsorbed O atoms, which we 
previously found to adsorb in non-perovskite sites^. 
This difference as well as the low mobility of paired 
hydroxyl groups could be important ingredients towards 
an atomic-scale understanding of the experimental 
reports that hydrogen and water increase the growth 
rate of SrTiOs.^i 
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